Whole-Genome Mapping of Epigenetic Modification of 5-Formylcytosine at Single-Base Resolution by Chemical Labeling Enrichment and Deamination Sequencing.

DNA cytosine methylation (5-methylcytosine, 5mC) is a predominant epigenetic modification that plays a critical role in a variety of biological and pathological processes in mammals. In active DNA demethylation, the 10-11 translocation (TET) dioxygenases can sequentially oxidize 5mC to generate three modified forms of cytosine, 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). Beyond being a demethylation intermediate, recent studies have shown that 5fC has regulatory functions in gene expression and chromatin organization. While some methods have been developed to detect 5fC, genome-wide mapping of 5fC at base resolution is still highly desirable. Herein, we propose a chemical labeling enrichment and deamination sequencing (CLED-seq) method for detecting 5fC in genomic DNA at single-base resolution. The CLED-seq method utilizes selective labeling and enrichment of 5fC-containing DNA fragments, followed by deamination mediated by apolipoprotein B mRNA-editing catalytic polypeptide-like 3A (APOBEC3A or A3A) and sequencing. In the CLED-seq process, while all C, 5mC, and 5hmC are interpreted as T during sequencing, 5fC is still read as C, enabling the precise detection of 5fC in DNA. Using the proposed CLED-seq method, we accomplished genome-wide mapping of 5fC in mouse embryonic stem cells. The mapping study revealed that promoter regions enriched with 5fC overlapped with H3K4me1, H3K4me3, and H3K27ac marks. These findings suggest a correlation between 5fC marks and active gene expression in mESCs. In conclusion, CLED-seq is a straightforward, bisulfite-free method that offers a valuable tool for detecting 5fC in genomes at a single-base resolution.

AlkB-Facilitated Demethylation Enables Quantitative and Site-Specific Detection of Dual Methylation of Adenosine in RNA.

RNA molecules undergo various chemical modifications that play critical roles in a wide range of biological processes. N6,N6-Dimethyladenosine (m6,6A) is a conserved RNA modification and is essential for the processing of rRNA. To gain a deeper understanding of the functions of m6,6A, site-specific and accurate quantification of this modification in RNA is indispensable. In this study, we developed an AlkB-facilitated demethylation (AD-m6,6A) method for the site-specific detection and quantification of m6,6A in RNA. The N6,N6-dimethyl groups in m6,6A can cause reverse transcription to stall at the m6,6A site, resulting in truncated cDNA. However, we found that Escherichia coli AlkB demethylase can effectively demethylate m6,6A in RNA, generating full-length cDNA from AlkB-treated RNA. By quantifying the amount of full-length cDNA produced using quantitative real-time PCR, we were able to achieve site-specific detection and quantification of m6,6A in RNA. Using the AD-m6,6A method, we successfully detected and quantified m6,6A at position 1851 of 18S rRNA and position 937 of mitochondrial 12S rRNA in human cells. Additionally, we found that the level of m6,6A at position 1007 of mitochondrial 12S rRNA was significantly reduced in lung tissues from sleep-deprived mice compared with control mice. Overall, the AD-m6,6A method provides a valuable tool for easy, accurate, quantitative, and site-specific detection of m6,6A in RNA, which can aid in uncovering the functions of m6,6A in human diseases.

Single-Nucleotide Resolution Mapping of N6-Methyladenine in Genomic DNA.

N6-Methyladenine (6mA) is a naturally occurring DNA modification in both prokaryotes and eukaryotes. Herein, we developed a deaminase-mediated sequencing (DM-seq) method for genome-wide mapping of 6mA at single-nucleotide resolution. The method capitalizes on the selective deamination of adenine, but not 6mA, in DNA mediated by an evolved adenine deaminase, ABE8e. By employing this method, we achieved genome-wide mapping of 6mA in Escherichia coli and in mammalian mitochondrial DNA (mtDNA) at single-nucleotide resolution. We found that the 6mA sites are mainly located in the GATC motif in the E. coli genome. We also identified 17 6mA sites in mtDNA of HepG2 cells, where all of the 6mA sites are distributed in the heavy strand of mtDNA. We envision that DM-seq will be a valuable tool for uncovering new functions of 6mA in DNA and for exploring its potential roles in mitochondria-related human diseases.

Single-Base Resolution Detection of N6-Methyladenosine in RNA by Adenosine Deamination Sequencing.

N6-Methyladenosine (m6A) is one of the most abundant and prevalent natural modifications occurring in diverse RNA species. m6A plays a wide range of roles in physiological and pathological processes. Revealing the functions of m6A relies on the faithful detection of individual m6A sites in RNA. However, developing a simple method for the single-base resolution detection of m6A is still a challenging task. Herein, we report an adenosine deamination sequencing (AD-seq) technique for the facile detection of m6A in RNA at single-base resolution. The AD-seq approach capitalizes on the selective deamination of adenosine, but not m6A, by the evolved tRNA adenosine deaminase (TadA) variant of TadA8e or the dimer protein of TadA-TadA8e. In AD-seq, adenosine is deaminated by TadA8e or TadA-TadA8e to form inosine, which pairs with cytidine and is read as guanosine in sequencing. m6A resists deamination due to the interference of the methyl group at the N6 position of adenosine. Thus, the m6A base pairs with thymine and is still read as adenosine in sequencing. The differential readouts from A and m6A in sequencing can achieve the single-base resolution detection of m6A in RNA. Application of the proposed AD-seq successfully identified individual m6A sites in Escherichia coli 23S rRNA. Taken together, the proposed AD-seq allows simple and cost-effective detection of m6A at single-base resolution in RNA, which provides a valuable tool to decipher the functions of m6A in RNA.

Enzymatic Cleavage-Mediated Extension Stalling Enables Accurate Recognition and Quantification of Locus-Specific Uracil Modification in DNA.

Chemical modifications in DNA have profound influences on the structures and functions of DNA. Uracil, a naturally occurring DNA modification, can originate from the deamination of cytosine or arise from misincorporation of dUTP into DNA during DNA replication. Uracil in DNA will imperil genomic stability due to their potential in producing detrimental mutations. An in-depth understanding of the functions of uracil modification requires the accurate determination of its site as well as content in genomes. Herein, we characterized that a new member of the uracil-DNA glycosylase (UDG) family enzyme (UdgX-H109S) could selectively cleave both uracil-containing single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA). Based on this unique property of UdgX-H109S, we developed an enzymatic cleavage-mediated extension stalling (ECES) method for the locus-specific detection and quantification of uracil in genomic DNA. In the ECES method, UdgX-H109S specifically recognizes and cleaves the N-glycosidic bond of uracil from dsDNA and generates an apurinic/apyrimidinic (AP) site, which could be broken by APE1 to form a one-nucleotide gap. The specific cleavage by UdgX-H109S is then evaluated and quantified by qPCR. With the developed ECES approach, we demonstrated that the level of uracil at position Chr4:50566961 in genomic DNA of breast cancer tissues was significantly decreased. Collectively, the ECES method has been proved to be accurate and reproducible in the locus-specific quantification of uracil in genomic DNA from biological and clinical samples.

Quantitative and site-specific detection of inosine modification in RNA by acrylonitrile labeling-mediated elongation stalling.

RNA molecules contain diverse modifications that play crucial roles in a wide variety of biological processes. Inosine is one of the most prevalent modifications in RNA and dysregulation of inosine is correlated with many human diseases. Herein, we established an acrylonitrile labeling-mediated elongation stalling (ALES) method for quantitative and site-specific detection of inosine in RNA from biological samples. In ALES method, inosine is selectively cyanoethylated with acrylonitrile to form N1-cyanoethylinosine (ce1I) through a Michael addition reaction. The N1-cyanoethyl group of ce1I compromises the hydrogen bond between ce1I and other nucleobases, leading to the stalling of reverse transcription at original inosine site. This specific property of stalling at inosine site could be evaluated by subsequent real-time quantitative PCR (qPCR). With the proposed ALES method, we found the significantly increased level of inosine at position Chr1:63117284 of Ino80dos RNA of multiple tissues from sleep-deprived mice compared to the control mice. This is the first report on the investigation of inosine modification in sleep-deprived mice, which may open up new direction for deciphering insomnia from RNA modifications. In addition, we found the decreased level of inosine at GluA2 Q/R site (Chr4:157336723) in glioma tissues, indicating the decreased level of inosine at GluA2 Q/R site may serve as potential indicator for the diagnosis of glioma. Taken together, the proposed ALES method is capable of quantitative and site-specific detection of inosine in RNA, which provides a valuable tool to uncover the functions of inosine in human diseases.

Engineered APOBEC3C Sequencing Enables Bisulfite-Free and Direct Detection of DNA Methylation at a Single-Base Resolution.

DNA methylation (5-methylcytosine, 5mC) is the most important epigenetic modification in mammals. Deciphering the roles of 5mC relies on the quantitative detection of 5mC at the single-base resolution. Bisulfite sequencing (BS-seq) is the most often employed technique for mapping 5mC in DNA. However, bisulfite treatment may cause serious degradation of input DNA due to the harsh reaction conditions. Here, we engineered the human apolipoprotein B mRNA-editing catalytic polypeptide-like 3C (A3C) protein to endow the engineered A3C (eA3C) protein with differential deamination activity toward cytosine and 5mC. By the virtue of the unique property of eA3C, we proposed an engineered A3C sequencing (EAC-seq) method for the bisulfite-free and quantitative mapping of 5mC in DNA at the single-base resolution. In EAC-seq, the eA3C protein can deaminate C but not 5mC, which is employed to differentiate C and 5mC in sequencing. Using the EAC-seq method, we quantitatively detected 5mC in genomic DNA of lung cancer tissue. In contrast to the harsh reaction conditions of BS-seq, which could lead to significant degradation of DNA, the whole procedure of EAC-seq is carried out under mild conditions, thereby preventing DNA damage. Taken together, the EAC-seq approach is bisulfite-free and straightforward, making it an invaluable tool for the quantitative detection of 5mC in limited DNA at the single-base resolution.

Whole-Genome Sequencing of 5-Hydroxymethylcytosine at Base Resolution by Bisulfite-Free Single-Step Deamination with Engineered Cytosine Deaminase.

The epigenetic modification 5-hydroxymethylcytosine (5hmC) plays a crucial role in the regulation of gene expression. Although some methods have been developed to detect 5hmC, direct genome-wide mapping of 5hmC at base resolution is still highly desirable. Herein, we proposed a single-step deamination sequencing (SSD-seq) method, designed to precisely map 5hmC across the genome at single-base resolution. SSD-seq takes advantage of a screened engineered human apolipoprotein B mRNA-editing catalytic polypeptide-like 3A (A3A) protein, known as eA3A-v10, to selectively deaminate cytosine (C) and 5-methylcytosine (5mC) but not 5hmC. During sequencing, the deaminated C and 5mC are converted to uracil (U) and thymine (T), read as T in the sequencing data. However, 5hmC remains unaffected by eA3A-v10 and is read as C during sequencing. Consequently, the presence of C in the sequence reads indicates the original 5hmC. We applied SSD-seq to generate a base-resolution map of 5hmC in human lung tissue. Our findings revealed that 5hmC was predominantly localized to CpG dinucleotides. Furthermore, the base-resolution map of 5hmC generated by SSD-seq demonstrated a strong correlation with prior ACE-seq results. The advantages of SSD-seq are its single-step process, absence of bisulfite treatment or DNA glycosylation, cost effectiveness, and ability to detect and quantify 5hmC directly at single-base resolution.

Bisulfite-Free and Single-Base Resolution Detection of Epigenetic DNA Modification of 5-Methylcytosine by Methyltransferase-Directed Labeling with APOBEC3A Deamination Sequencing.

DNA methylation (5-methylcytosine, 5mC) is the most prevalent epigenetic modification that is predominantly found in CG dinucleotides in mammalian genomes. In-depth investigation of the functions of 5mC heavily relies on the quantitative measurement of 5mC at single-base resolution in genomes. Here, we proposed a methyltransferase-directed labeling with APOBEC3A (A3A) deamination sequencing (MLAD-seq) method for the single-base resolution and quantitative detection of 5mC in DNA. In MLAD-seq, a mutant of DNA methyltransferase, M.MpeI-N374K, is utilized to selectively transfer a carboxymethyl group to the 5 position of cytosine in the CG dinucleotide to form 5-carboxymethylcytosine (5camC) using carboxy-S-adenosyl-l-methionine (caSAM) as the cofactor. After A3A treatment, 5camC is resistant to the deamination and base pairs with guanine. Thus, the cytosines in CG sites are read as C in sequencing. On the contrary, the methyl group in 5mC inhibits its carboxymethylcytosine by M.MpeI-N374K and therefore is readily deaminated by A3A to produce thymine that pairs with adenine and is read as T in sequencing. The differential readouts from C and 5mC in the MLAD-seq enable the single-base resolution mapping of 5mC in CG sites in DNA. With the developed MLAD-seq method, we observed the hypermethylation in the promoter region of retinoic acid receptor ß (RARB) gene from human nonsmall cell lung tumor tissue. Compared to harsh reaction conditions in bisulfite sequencing that could lead to significant degradation of DNA, the whole procedure of MLAD-seq is carried out under mild conditions, which will avoid DNA damage. Thus, MLAD-seq is more suitable in the scenario where only limited input DNA is available. Taken together, the MLAD-seq offers a valuable tool for bisulfite-free, single-base resolution and quantitative detection of 5mC in limited DNA.

Genome-wide mapping of N 4-methylcytosine at single-base resolution by APOBEC3A-mediated deamination sequencing.

N 4-methylcytosine (4mC) is a natural DNA modification occurring in thermophiles and plays important roles in restriction-modification (R-M) systems in bacterial genomes. However, the precise location and sequence context of 4mC in the whole genome are limited. In this study, we developed an APOBEC3A-mediated deamination sequencing (4mC-AMD-seq) method for genome-wide mapping of 4mC at single-base resolution. In the 4mC-AMD-seq method, cytosine and 5-methylcytosine (5mC) are deaminated by APOBEC3A (A3A) protein to generate uracil and thymine, both of which are read as thymine in sequencing, while 4mC is resistant to deamination and therefore read as cytosine. Thus, the readouts of cytosines from sequencing could manifest the original 4mC sites in genomes. With the 4mC-AMD-seq method, we achieved the genome-wide mapping of 4mC in Deinococcus radiodurans (D. radiodurans). In addition, we confirmed that 4mC, but not 5mC, was the major modification in the D. radiodurans genome. We identified 1586 4mC sites in the genome of D. radiodurans, among which 564 sites were located in the CCGCGG motif. The average methylation levels in the CCGCGG motif and non-CCGCGG sequence were 70.0% and 22.8%, respectively. We envision that the 4mC-AMD-seq method will facilitate the investigation of 4mC functions, including the 4mC-involved R-M systems, in uncharacterized but potentially useful strains.

Bisulfite-free and single-nucleotide resolution sequencing of DNA epigenetic modification of 5-hydroxymethylcytosine using engineered deaminase.

The discovery of 5-hydroxymethylcytosine (5hmC) in mammalian genomes is a landmark in epigenomics study. Similar to 5-methylcytosine (5mC), 5hmC is viewed as a critical epigenetic modification. Deciphering the functions of 5hmC necessitates the location analysis of 5hmC in genomes. Here, we proposed an engineered deaminase-mediated sequencing (EDM-seq) method for the quantitative detection of 5hmC in DNA at single-nucleotide resolution. This method capitalizes on the engineered human apolipoprotein B mRNA-editing catalytic polypeptide-like 3A (A3A) protein to produce differential deamination activity toward cytosine, 5mC, and 5hmC. In EDM-seq, the engineered A3A (eA3A) protein can deaminate C and 5mC but not 5hmC. The original C and 5mC in DNA are deaminated by eA3A to form U and T, both of which are read as T during sequencing, while 5hmC is resistant to deamination by eA3A and is still read as C during sequencing. Therefore, the remaining C in the sequence manifests the original 5hmC. By EDM-seq, we achieved the quantitative detection of 5hmC in genomic DNA of lung cancer tissue. The EDM-seq method is bisulfite-free and does not require DNA glycosylation or chemical treatment, which offers a valuable tool for the straightforward and quantitative detection of 5hmC in DNA at single-nucleotide resolution.

Single-Base Resolution Detection of Adenosine-to-Inosine RNA Editing by Endonuclease-Mediated Sequencing.

RNA molecules contain diverse modifications that play crucial roles in a wide variety of biological processes. Adenosine-to-inosine (A-to-Ino) RNA editing is one of the most prevalent modifications among all types of RNA. Abnormal A-to-InoRNA editing has been demonstrated to be associated with many human diseases. Identification of A-to-Ino editing sites is indispensable to deciphering their biological roles. Herein, by employing the unique property of human endonuclease V (hEndoV), we proposed a hEndoV-mediated sequencing (hEndoV-seq) method for the single-base resolution detection of A-to-InoRNA editing sites. In this approach, the terminal 3'OH of RNA is first blocked by 3'-deoxyadenosine (3'-deoxy-A). Specific cleavage of Ino sites by hEndoV protein produces new terminal 3'OH, which can be identified by sequencing analysis, and therefore offers the site-specific detection of Ino in RNA. The principle of hEndoV-seq is straightforward and the analytical procedure is simple. No chemical reaction is involved in the sequencing library preparation. The whole procedure in hEndoV-seq is carried out under mild conditions and RNA is not prone to degradation. Taken together, the proposed hEndoV-seq method is capable of site-specific identification of A-to-Ino editing in RNA, which provides a valuable tool for elucidating the functions of A-to-Ino editing in RNA.

Direct decarboxylation of ten-eleven translocation-produced 5-carboxylcytosine in mammalian genomes forms a new mechanism for active DNA demethylation.

DNA cytosine methylation (5-methylcytosine, 5mC) is the most important epigenetic mark in higher eukaryotes. 5mC in genomes is dynamically controlled by writers and erasers. DNA (cytosine-5)-methyltransferases (DNMTs) are responsible for the generation and maintenance of 5mC in genomes. Active demethylation of 5-methylcytosine (5mC) is achieved by ten-eleven translocation (TET) dioxygenase-mediated oxidation of 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). 5fC and 5caC are further processed by thymine DNA glycosylase (TDG)-initiated base excision repair (BER) to restore unmodified cytosines. The TET-TDG-BER pathway could cause the production of DNA strand breaks and therefore jeopardize the integrity of genomes. Here, we investigated the direct decarboxylation of 5caC in mammalian genomes by using metabolic labeling with 2'-fluorinated 5caC (F-5caC) and mass spectrometry analysis. Our results clearly demonstrated the decarboxylation of 5caC occurring in mammalian genomes, which unveiled that, in addition to the TET-TDG-BER pathway, the direct decarboxylation of TET-produced 5caC constituted a new pathway for active demethylation of 5mC in mammalian genomes.

Single-Nucleotide Resolution Analysis of 5-Hydroxymethylcytosine in DNA by Enzyme-Mediated Deamination in Combination with Sequencing.

The report of the existence of 5-hydroxymethylcytosine (hm5C) in mammalian genomes is a milestone discovery. hm5C is now generally viewed as the sixth base of DNA with important functions on epigenetic regulation. The in-depth investigation of the biological functions of hm5C requires elucidating the distribution patterns of hm5C in genomes, better in single-nucleotide resolution. It was reported that the cytosine deaminases of the APOBEC (apolipoprotein B mRNA-editing catalytic polypeptide-like) family are nucleic acid editing enzymes and can deaminate cytosine (C) to form uracil (U). Particularly, a subfamily of APOBEC (APOBEC3A) can efficiently deaminate both C and 5-methylcytosine (m5C). In the current study, we identified that APOBEC3A protein can effectively deaminate C, m5C, and hm5C but shows no observable deamination activity toward glycosylated hm5C (ß-glucosyl-5-hydroxymethyl-2'-deoxycytidine, ghm5C) by using the restriction enzyme-based assay and liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) analysis. By virtue of the differential deamination activity of APOBEC3A toward C, m5C, and ghm5C in conjugation with sequencing, we developed the single-nucleotide resolution analysis of hm5C in DNA. In this analytical strategy, the original C and m5C in DNA will be deaminated by APOBEC3A to form U and thymine (T), both of which will read as T during sequencing, while ghm5C is resistant to deamination and will read as C during sequencing. Therefore, the remaining C in the sequence context only could come from original hm5C, which offers the single-nucleotide resolution analysis of hm5C in DNA. This APOBEC3A-mediated deamination sequencing (AMD-seq) is straightforward and involves no bisulfite treatment, which avoids the substantial degradation of DNA. Future application of this strategy can be performed for the reliable mapping of hm5C in genome-wide scale at the single-nucleotide resolution.